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ABSTRACT 

We have used the ACS camera on HST to obtain (V, I) photometry for 57,000 
red-giant stars in the halo of the Leo elliptical NGC 3377, an intermediate- 
luminosity elliptical. We use this sample of stars to derive the metallicity distri- 
bution function (MDF) for its halo field stars, and comment on its chemical evo- 
lution history compared with both larger and smaller E galaxies. Our ACS /WFC 
field spans a radial range extending from 4 to 18 kpc projected distance from the 
center of NGC 3377 and thus covers a significant portion of this galaxy's halo. 
We find that the MDF is broad, reaching a peak at log (Z/Z Q ) ~ —0.6, but 
containing virtually no stars more metal-poor than log (Z/Z Q ) = —1.5. It may, 
in addition, have relatively few stars more metal-rich than log (Z/Zq) = —0.3, 
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although interpretation of the high-metallicity end of the MDF is limited by pho- 
tometric completeness that affects the detection of the reddest, most metal-rich 
stars. NGC 3377 appears to have an enrichment history intermediate between 
those of normal dwarf ellipticals and the much larger giants. As yet, we find no 
clear evidence that the halo of NGC 3377 contains a significant population of 
"young" (< 3 Gy) stars. 

Subject headings: galaxies: elliptical — galaxies: individual (NGC 3377) 



Introduction 



The imaging cameras on the HST have provided powerful tools to study large pho- 
tometric samples of individual stars in nearby galaxies that were quite out of reach from 
ground-based instruments. Although even the best photometry cannot provide the level of 
detail that is achievable by high-resolution spectroscopy, these deep stellar samples open a 
valuable route to understanding the early chemical evolution of their host galaxies, compared 
with the much cruder information from integrated-light photometry and spectroscopy. 

Among the possible target galaxies within reach of the HST, the ellipticals are of spe- 
cial interest because they may result from the widest possible range of formation histories, 
from hierarchical merging at very early times, to recent major mergers, to later growth 
by satellite accretion. Representative dwarf ellipticals are readily accessible within the Lo- 
cal Group (particularly NGC 147, 185, 205), but the Local Group gives us no examples 
of giant or intermediate-sized E galaxies to work with. In a series of previo us papers 
flHarris. Harris. & Poole] Il999l : lHarris fc Harrisl boool . I2OO2I : iReikuba et alTl2005h . we pre- 
sented analyses of samples of stars in the closest giant elliptical, NGC 5128. These studies 
cover four locations in its halo at projected distances ranging from 8 to 40 kp c. In all four of 



these studies the red-giant stars are clearly resolved, and in the deepest one (IRejkuba et al. 



20051 ) . even the old horizontal-branch population is reached. We found that in all four lo- 
cations in the halo of this giant elliptical, the metallicity distribution of the giant stars is 



1 Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. These observations are associated with program #9811. Support 
for this work was provided in part by NASA through grant number HST-GO-09811.01-A from the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS 5-26555. 
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clearly metal-rich (with ( Fe/H ) ~ —0.4) and broad, but with extremely small proportions 
of classically metal-poor stars (i.e. those with [Fe/H] < — 1). 

Generalizing the results for NGC 5128 to all large ellipticals may, however, be com- 
promised by the possibility that NGC 5128 could be a major merger remnant in which the 
majority of its halo stars were actually formed within the metal-rich disks of its progenitor 
galaxies. In fact, both a major-merger or hierarchical-merging approach are ca pable of cre- 



ating an MDF in the ha lo with the same basic characteristics as we observe ( iBekki et al. 



20031 : iBeasley et al.ll2003l ). To extend the range of information we have to work with, clearly 
it is of great interest to obtain the metallicity distribution function (MDF), based directly 
on samples of individual stars, in other galaxies of the widest possible range of properties. 

NGC 5128, at a distance of only 4 Mpc, is by far the closest easily accessible big elliptical. 
To reach other ellipticals, we must step to the Leo group objects at D ~ 10 Mpc, and then 
outward to Virgo (D ~ 16 Mpc) and Fornax (D ~ 19 Mpc), which have the nearest accessible 
large samples of target galaxies. The technological gains afforded by the step from the HST 
WFPC2 camera to the newer ACS, with its much higher sensitivity, spatial resolution, and 
field of view, now bri ng these more d i stant targets within reach (see, for example, the recent 



photometric work of IWilliams et al.l (120071 ) to resolve the brightest ~ 1.5 mag of the red- 
giant branches in four globular clusters at the Virgo distance). Roughly speaking, observing 
the halo of NGC 5128 with WFPC2 is comparable (for the same number of HST orbits) 
with observing the Leo group ellipticals with ACS. 

In this paper, we present new color-magnitude photometry for the halo stars in the 
intermediate-sized E5 elliptical NGC 3377 in the Leo group, and briefly compare its metal- 
licity distribution function (MDF) with those of both a dwarf and a giant. NGC 3377 has 
a luminosity of My ~ —19.9 falling between those of the Local Group dwarfs and the giant 
NGC 5128 that have been observed in previous studies; it therefore provides an attractive 
bridge in the parameter space of galaxy properties that we can add to the discussion of 
chemical evolution models. The key parameters for the galaxy are summarized in Table 1. 



2. Observations and Data Reduction 

The imaging data we use for this study were obtained in our HST program 9811. NGC 
3377 is at relatively high galactic latitude and low foreground reddening {b = 58?3, E B _ V 
= 0.03) and is thus well placed for deep photometry of its halo stars, as free from field 
contamination and foreground absorption as we can hope to have. 

We used the Advanced Camera for Surveys in its Wide Field Channel, with image scale 
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0'/05 per pixel, to reach the maximum photometric depth and field coverage. A single field 
was targeted at a = 10 A 47 m 49 m 00, 5 = + 13?55^40"0 (J2000). This field is 3!82 southeast of 
the galaxy center, equivalent to ~12 kpc projected linear radius. Since the effective radius 
of the galaxy light profile is R e = l'.l, our field is well beyond the central bulge and can 
plausibly be considered as giving us a first look at the genuine halo. Fields further in would 
obviously be of interest to track any metallicity gradient that the system might have, but 
would also be much more difficult to measure because of considerably increased crowding. 

The filters used were the ACS/WFC "wide V" (F60QW) and "wide I" (F814W), the 
same ones as in our previous NGC 5128 studies. The (V — I) color index is a very useful 
metallicity indicator particularly for old red-giant stars in the range [Fe/H] ~ — 1.5 to ~ 0, 
which as will be seen below, is our main range of interest. Total exposure times were 38500 
sec over 15 full-orbit exposures in FQOQW, and 22260 sec over 9 full-orbit exposures in 
F81AW, adding up to a total of 24 orbits on this single field. The individual exposures were 
dithered over steps from to 20 pixels, allowing elimination of most cosmic rays, bad pixels, 
and other artifacts on the detector. To prepare the images for photometry, we extracted 
the drizzled individual images from the HST Data Archive, registered them to within 0.05 
pixels, and median-combined them. This procedure gave us a single very deep exposure in 
each filter. In Figure 1, we show the combined J— band image, demonstrating that the halo 
is well resolved into stars throughout. A small inset portion of the field is shown in Figure 
2. Many faint background galaxies can be seen through the halo, but the overwhelming 
majority of detected objects in this field are the red giant stars of NGC 3377 itself. 

To carry out the photometry on the two combined frames we used the standalone version 
of DAOPHOT in its latest (daophot 4) version as written by one of us (PBS). The normal 
sequence of find/phot /allstar was used, with point spread functions defined from the average 
of 50 to 60 stars spread across the field. The pixel coordinate systems on the two images were 
registered, and then the detected stars on each one were matched up to find those measured 
in both colors. The number density of stars on the frame was reasonably high (the average 
separation between detected stars is ~ 7 px, while the FWHM of the point spread function 
is 2.3 px or 0'.'115); also, differences in magnitude are not useful for deciding on matchup 
of stars between the V and / frames since the target stars have a very large true range in 
color. Matchup was therefore done by coordinates alone in a careful iterative procedure: five 
iterations were done with successively increasing matching radius, starting with Ar = 0.2 
px and increasing to 1.5 px. In practice, more than 95% of the 75000 matches between the 
V and I images fell well within 1 pixel of each other. 

The next step was to define the cleanest possible subset of these matches. We rejected 
any detected objects that fell in the regions of 5 bigger-than-average background galaxies 
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and one bright star, simply by masking out those regions. We also rejected objects with 
very poor goodness-of-fit to the PSF (xv > 2.5, xi > 2.0); and any with much larger than 
normal ALLSTAR magnitude uncertainties (ey, ej > 0.25). 

Although there is a very noticeable gradient in the density of stars from the lower right 
portion of the field to the upper left (see Fig.l), none of the areas are excessively crowded 
to the point where photometry via normal PSF fitting is unusually difficult. Our data do 
not enter the high-crowding regime that can be found in other examples in the literature, 
such as (among others) the i nner halo of NGC 51 28 ([Harris fc Harris! |2002| ). the bulge of 
the Local Group d warf M32 (iGrillmair et al.lll996h . or the halo of the Virgo dE VCC1104 
(IHarris et al.lll998l ). Nevertheless, to extract a conservatively "uncrowded" sample of stars 
to work with, we calculated the nearest-neighbor distance (NND) for every measured star on 
both frames. We define the NND for a star simply as the distance to the closest neighboring 
star on the list regardless of relative brightness. The distribution of the NND values is shown 
in Figure [3j showing that the typical separation between objects is ~ 5 px or about twice 
the FWHM. To cull the photometry files of stars that are even mildly crowded, we decided 
to reject any object for which NND < 3 px, on either the V or I frames. This step removed 
about 17% of the I detections and 12% of the V detections. We regard the NND cutoff of 3 
px as an extremely stringent choice; in fact, allstar has no difficulty handling cases down to 
NND ~ FWHM or even less. 

The final result of the photometry, after removal of (a) the masked regions, (b) objects 
with high x an d high magnitude uncertainties, and (c) stars for which any concern existed 
about crowding, left a total of 57039 well measured, uncrowded stars. 

For calibration of the photometry, we chose to transform the filter-bas ed magnitudes 

F606W, F814W into V, I in order to facilitate comparis ons with previous work (IHarris. Harris, fc Poole 
19991 ; IHarris fc Harris! |2000| . 120021 ; iRejkuba et al.ll2005l ). First, we used the daophot /substar 
routine to subtract all but the brighter stars from the final averaged pair of images. We 
then obtained aperture photometry of these now-isolated bright stars to correct the allstar 
instrumental magnitudes to aperture mag nitudes, and extrapo lated these to large radius 
following the standardized prescriptions of ISirianni et al. These gave us the mag- 

nitudes FQ0QW and F8UW on the natural ACS VEGAMAG filter system. Finally, we 
independently derived transformations of these to the standard VI system by measuring an 
extensive set of images of the NGC 2419 standard field from the HST Archive, and compar- 
ing these with ground-based stan dard data in t he same field. The ground-based data are 
part of PBS's ongoing programme (jStetson!l2000l ) t o maintain and upgrade an all-sky system 



of photometric standards on the BVRI sy stem of 



NGC 2419 field in particular is discussed by lStetsonl (120051 ) . This remote globular cluster has 



Landqltl (119921 ) . The photometry of the 
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the advantages (for calibration purposes) of a blue horizontal branch at V~21, a red giant 
branch tip at (V, V-I) = (17.3, 1.50), and a wealth of foreground stars. The 1,257 stars 
adopted as standards in this field have a median of 46 groundbased observations per star in 
each of the V and I filters, a median standard error of 0.0054 mag in V and 0.0058 mag in 
J, and a magnitude and color range 12.66 < V < 22.64, —0.05 < (V — I) < 3.04. Keeping 
quadratic terms in color index, we found for these 
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These equations reproduce the NGC 2419 standard stars in all bands to within a scatter 
a = ±0.03 mag. The equation for B is not used in this pap er but is listed for in formation. 
These transformations are very close to those publ ished in ISirianni et al.l (120051 ). and also 
to those derived for our NGC 5128 outer-halo field (IRejkuba et al.ll2005l ). which were taken 
during the same Cycle. The final color-magnitude diagram for the complete sample of 57039 
stars is shown in Figure H] and will be discussed below. 

We estimated the internal photometric uncertainties and the detection completeness 
through a series of artificial-star tests with the addstar component of DAOPHOT. Stars 
were added to the combined V and / images in groups of 1000 over a wide range of mag- 
nitudes; these experiments were done independently on the F606W and F814W images. 
The images were remeasured in the same way as the original frames. The fraction / of 
stars recovered, as a function of instrumental magnitude, is shown in Figure [5J (The mag- 
nitudes here are the filter-based ones F606W and F814W, discussed below.) The limits 
of our data, defined as the magnitudes at which / drops to 0.5, are F606W(lim) = 28.95 
and F8l4W(lim) = 27.70. The tren d of / with magnitu de is well described in each case 
by a Pritchet interpolation function ([Fleming et al.lll995l ). an analytic function with two 
parameters: the limiting magnitude, and a parameter a giving the steepness of the dropoff. 
For these images we find a = 2.5 for V and a = 2.7 for J. Because of the color terms in 
the transformations (Eqs. 2,3), these limits in the native filter-defined magnitudes do not 
correspond to single V or I values, but for a giant star with a typical color of (V — I) — 2, 
the limits are V(lim) ~ 29.4 and I(lim) ~ 27.7. 

In the color-magnitude diagram of Fig. HJ the 50% completeness lines are shown. An 
important feature of these lines is that the limiting curve for V (the upward- slanting line 
on the right side of the CMD) cuts off our ability to see any extremely red stars that might 
actually be present; these would fall at the most metal- rich end of our metallicity distribution 
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function. Considerably deeper exposures in V will be needed to explore the true "red limit" 
of the giant stars in this galaxy. Within the limits imposed by the photometry, we explicitly 
take into account the completeness fraction / in our derivation (below) of the metallicity 
distribution. As will be seen later, the completeness cutoff may affect how much of the 
metal-rich end of the MDF is ultimately detectable. 

We find that the 50% completeness level does not change significantly with radius R 
from galaxy center, except perhaps marginally in the region R < 2' of highest crowding. 
A visual confirmation is shown in Figure El where we subdivide the CMD data into four 
radial regions. The placements of the completeness cutoffs do not change relative to the 
distribution of stars (we will further justify this point in the later discussion on the radial 
gradient of the metallicity distribution). A further indication of the relative importance of 
crowding is shown in Figure where we show the local sky intensity around each measured 
star (in digital units, directly as returned by allstar) as a function of radius from the center 
of NGC 3377. The inward increase becomes steeper within R = 2', but the net change over 
the entire frame is quite modest. However, in the discussion that follows (for the distance 
measurement and the metallicity distribution function), to remove any residual concerns 
about crowding issues we further restrict our analysis to the "safest" region R > 2'. 

The random uncertainties of the photometry, as derived from the artificial-star tests, 
are shown in Figure [HJ The rms uncertainty rises smoothly with magnitude in a roughly 
exponential manner, shown by the curves in the Figure. These have equations cr(V) = 0.05 
exp((F606H/ - 26.00)/1.8) and er(J) = 0.05 exp((F814W - 25.5)/1.3). In the analysis, we 
make no use of the data fainter than the 50% completeness level (marked by the vertical 
line in the figures). Finally, in Figure [9] we show the mean trends for systematic bias in the 
photometry: a sample of the artificial-star tests is plotted, showing the median difference 
Am = (input — measured) for both filters. The median lines, plotted in 0.1-mag bins, 
show a consistent trend for stars to be measured 0.01 — 0.02 mag too bright, a common 
feature of photometry of faint stars within moderately crowded fields. Fainter than the 50% 
completeness level, the mean curves diverge strongly to positive Am and the photometry 
becomes systematically unreliable. It is important to note, however, that for the stars 
brighter than the completeness limit, the mean biasses in the two filters run closely parallel, 
and the net bias in mean color (V — I) - an important consideration since the color index 
is our main metallicity indicator - is negligibly small and has no effect on the derivation of 
the MDF. In addition, the random uncertainties have no significant effect on the inferred 
properties of the MDF, because we use only the brightest ~ 1 mag of the RGB where the 
intrinsic color range of the stars is far larger than the photometric scatter. At the very top 
of the RGB, the intrinsic spread is A (V — I) ~ 1.6 mag while the photometric scatter is 
±0.1 mag or less. 
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Distance Calibration 



NGC 3377 is part of the Leo I group along with several other large galaxies, most of 
which are spiral-type. Numerous measurements of distance to these individual Leo members 
are in the literature from a variety of well established distance indicators including Cepheids, 
planetary nebula luminosity function (PNLF), surface brightness fluctuation (SBF), and the 
tip of the old- red-giant branch (TRGB). 



NGC 3377: For NGC 3377 itself, t he SBF method (ITonrv et al.l l200lh gives (m - M) 
H = 30.25, while the PNLF method Jciardullo et alill989h gives^ = 30.07 ± 0.18. 



NGC 3379: For the single giant elliptical in Leo I, the TRGB method has been applied 
through HST pho tometry in both the optical I band and t he near-infrared. The former 
JSakai et all Il997l ) gives \i = 30 30 ± 0.27, and the l atter (iGregg et all l2004h gives y, = 
30.17±0 12. The PNLF m ethod Jciardullo et al.lll989h gives // = 29.96 ±0.16, and the SBF 
method jTonrv et alJl200lh // = 30.12. 



Large Spirals: Cepheid-based distances have been publish ed for some of the m ajor spirals in 
the group. The se include u = 30.0 1 ± 0.19 for NGC 3351 jGraham et alJll997h . 30.25 ±0.18 
for NGC 3368 ([Tanvir et aDll999h . and 30.10 ± 14 (un weighted average of 27 Cepheids; 
Saha et alill999h or 29.71 ± 0.08 (IFreedman et al.l boOlh for NGC 3627. PNLF distances 
for these galaxies incl ude y = 30.05 ± 0.16 for NGC 3351 and 29.99 ± 0.08 for NGC 3627 
Jciardullo et alil2002h . 



Treating all these Leo members as if they are at the same true distance from us and taking 
an indiscriminate average of all these measurements suggests a rough consensus near \i ~ 
30.1 ± 0.05, or D ~ 10.4 Mpc for the group as a whole. The galaxy-to-galaxy dispersion 
of these measurements, which is = 0.17 mag, is quite similar to the typical internal 
uncertainties of each one and gives no strong evidence that the distance depth of the group 
is an important factor. 

Our new photometry penetrates well into the the old-halo red giant branch of NGC 
3377, with a large sample of stars, and thus provides a new opportunity to use the TRGB 
distance indicator rather precisely. The key parameter is the "tip magnitude" which is, 
physically, the luminosity of the helium flash in the core of the red giant as it reaches the 
top of its first ascent along the giant branch. The bolometric luminosity of the RGB tip is, 
fortunately, only mildly dependent on metallicity for old stars, allowing it to be turned into 
an accurate standard candle. 



Our analysis follows the methods used in lSakai et al.l (119961 . 119971 ) and lHarris. Harris, fc Poole 
(119991 ) . among others: we plot up the luminosity function of the RGB stars and use the de- 
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tailed shape of the LF in the / band to define the onset of the RGB. For stars more metal-poor 
than [Fe/H] ~ —0.7 (which include the majority of the ones we measure here; see next sec- 
tion), the I band has the strong advantage that the differential bolometric correction across 
the top of the RGB is almost cancelled by the dependence of M^itip) on metallicity, leaving 
Mj(tip) with only a gradual slope with increasing color. 

The luminosity function is shown in Figure [10], based on 49380 stars beyond R > 2' 
from the galaxy center. It has been smoothed with a Gaussian kernel of crj = 0.02 mag. 
Completeness corrections are quite unimportant here, since the / = 0.5 completeness level 
is I ~ 27.6, much fainter than the well resolved top of the RGB. The fact that the LF rises 
gradually upward at the TRGB (rather than the ideal case of an abrupt jump) is a normal 
consequence of an intrinsically steep LF at that point, convolved with the observational 
photometric scatter, and adding in the smoothing effect of the slight downward slope of the 
tip toward redder colors, and a small a mount of field contamination, all of which act to b lur 



out the TRGB point to some extent (see lHarris et al.lll998l ; lHarris. Harris, fc Pooldll999l . for 



detailed discussion and methodology). The actual tip is clearly somewhere near I = 26, but 
to define the true TRGB precisely, we look for a strong change in the slope of the LF as we 
go to fainter magnitudes across the tip. The numerically calculated first derivative of the LF, 
dn/dl, is shown in the second panel. A sharp upturn is present past I = 26 with well defined 
peaks that appear at / = 26.2 and 26.3. These are independent of the precise value of the 
smoothing kernel within broad limits. To decide which of these to pick we use the "edge 
response filter" or ERF (the numerical second derivative of the LF), plotted in the lower 
panel of Fig. [TDJ Here, we adopt the first and most prominent peak at I (tip) = 26.2 ± 0.1. 

The distance modulus follows immediately once we apply a fiducial value for Mi(tip). 
The most well establishe d recent calibration fr om within the Milky Way is from a very large 



sample of stars in u Cen (IBellazzini et al.ll2004l). giving M jitiv) = —4.05 ±0.12. Calibrations 



based on theoretical RGB models (e.g. ISalaris et al.ll2002l ) yield values in the range —3.95 to 



—4.22 depending on both the details of the stellar physics and on observational constraints, 
but are entirely consistent with the u Cen value, which we adopt here. We therefore obtain 
(m — M)i = 30.25 ±0.15 for NGC 3377. This must be corrected for a foreground absorption 
of Aj = 0.07 ± 0.02, giving a final TRGB distance measurement \x = 30.18 ± 0.16. 

Averaging the TRGB distance in with the SBF and PNLF measurements listed above, 
and giving the three methods equal weights, we arrive at an average (m — M)q = 30.17±0.10, 
or D = 10.8 ±0.5 Mpc for NGC 3377 itself. These three methods - all based, necessarily, 
on the properties of the old stellar population - are in excellent mutual agreement, and we 
consider the distance to this system to be as well established as any in the local galactic 
neighborhood. 
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4. The Nature of the Brighter Stars 



In Figs. H] and [6] there are numerous stars scattered above the TRGB. There are Nb rig ht — 
1360 such stars within 24 < / < 26. Could these be a signature of a younger, intermediate- 
age population (r < 3 Gy), which would indicate some more recent star formation in the 
galaxy's history? in Figure [TTJ we show the spatial distribution of these stars across our 
ACS field, compared with the distribution of the stars in the brightest half-magnitude of the 
RGB itself. Both types of stars show an obvious gradient decreasing from galaxy center, and 
so a significant fraction of Nbright must be genuine members of NGC 3377. 

Before placing limits on the number of such stars that might genuinely be present, we 
need to rule out other possible contributors, including (a) field contamination, (b) accidental 
blends of two normal RGB stars, or (c) giants or AGB stars in highly evolved or temporary 
states including long-period variables (LPVs), w hich can contribute noticeable number s of 
supra-TRGB stars even in an old population (see lRejkuba et al.ll2003l ; iGregg et al.ll2004 , for 
more exten sive discussion and modelling) . We dealt with a similar issue for the halo stars in 
NGC 5128 ( Harris. Harris, fc Pooldll999l ) and found that only ~1% of the population there 
could be ascribed to the intermediate-age category. 

The first option (field-contaminating objects) will include both foreground stars and 
faint, very small background galaxies that are near-starlike in appearance. The number 
of foreground stars in the direction of NGC 3377 and over the area of ACS/WFC should 



be N < 20 from Galactic models (e.g. iBahcall &: Soneiral Il98ll ). but the number of faint 



misidentified background galaxies is almost certainly larger and would best be measured 
from an adjacent 'control' field, which we do not have. However, we can specify a reasonable 
upper limit to Nf idd by using the outermost regions of our measured field. For R GC > 4f5 
(an area of 2.78 arcmin 2 ) there are 94 objects within 24 < / < 26, suggesting that over the 
entire 11.33 arcmin 2 field Nf ie id should be at most 376. The number density of stars is still 
declining at the outer edges of the field, and so this number must be a generous upper limit. 
Since we see 1360 such objects over the whole field, this suggests ~ 1000 — 1200 of them are 
intrinsic to NGC 3377. 

Some objects could appear in this bright range because of accidental blends of two stars 
that are both near the upper end of the RGB. Blends like this are very likely to be responsible 
for much of the scatter of stars appearing just above the nominal RGB tipE, but not for the 



2 There are ~ 1300 stars in our CMD in the small interval / = 24.0 - 24.2 just above the RGB tip that 
are responsible for blurring out the definition of the TRGB and creating the smooth rolloff in the luminosity 
function just above the tip; see Figure 10. Most of these objects can be understood as due to blends of the 
brighter RGB stars with the huge number of RGB stars > 1.5 mag fainter than the tip. 
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ones ~ 0.7 mag or more above it. Statistically the number of blended pairs will increase as 
the square of the number of stars per unit area, 



N, 



blend 



2 d 2 



(4) 



where N+ is the number of stars on the frame capable of generating a blended pair brighter 
than the TRGB, q is the radius of one resolution element, and d 2 is the area of the field. 
Adopting q ~ 2 px, and using the fact that there are ~ 30000 stars within 1 magnitude 
of the RGB tip, we expect N^end ~ 350, or about one-quarter of all the objects appearing 
above the TRGB. 



sec 



The presence of LPV -type stars is favored in an intrinsically metal-rich old population 



Rejkuba et al.ll2003l ). and their numbers should be proportional to the total luminosity 



of the whole population contained in our field. Extrapolating from the stars brighter than 
our completeness limit with a standard Population II luminosity function, we estimate very 
roughly that V(int) = 14.7 or L ~ 8 x 10 7 L o . From|Renzinj (119981 ) we then find N LPV ~ 400, 
similar to the expected number of accidental blends. An observational confirmation that 



LPV-type variables should indeed be there to be found is discussed by iGregg et al.l (12004 ). 
who present HST/NICMOS near-infrared photometry of the RGB stars in inner-halo fields 
for NGC 3379, the other Leo elliptical. They find (see their Table 3) that between one-quarter 
and one-h alf of the detec t ed sta rs above the TRGB are variable, consistent with our rougher 



estimate. 



Rejkuba et al.l (120031 ) found a total of more than 1100 LPVs in two fields around 
NGC 5128 covering an area of 10.46 arcmin 2 , and state that these make up roughly half of 
the total number of stars above the TRGB, again consistent with the estimates made here. 

In summary, we suggest that more than half the population of stars that lie clearly 
above the RGB tip are due to a combination of blends, temporary high-luminosity states 
of the giants, and a small amount of field contamination. The remainder not accounted 
for by these effects is thus N bright ~ 400 ± 100, or ~1% of our observed total on the field. 
We conclude tentatively that there is no evidence for any significant "young" (r < 3 Gy) 
population of stars in our field. 

These rough estimates are not intended to replace a more comprehensive population 
synthesis analysis, based on simulating the observed RGB with sta rs drawn from model 
isochrones over a complete range of ages and metallicities (see, e.g. 



Williams et al. 2007 



for a recent example). In a later stage of this work we will develop such simulations and 
constrain the age range of the halo stars more completely. 
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5. The Metallicity Distribution 



5.1. Deriving the MDF 

We are now in a position to use the colors of the RGB stars to gain a first direct 
measurement of the metallicity distribution function (MDF) in this galaxy. To enable direct 
comparisons with other syste ms, we follow the same method as used previously for NGC 5128 
and selected dwarf ellipticals (lHarris. Harris, fc Pooldll999l : lHarris fc Harris! l2000l . 120021 ). A 
finely spaced grid of RGB evolutionary tracks is superimposed on the measured CMD, and 
we interpolate between these to tag every star with a heavy-element abu ndance Z based on 



i ts loc ation within the tracks. As before, we use the a— enhanced tracks of IVandenBerg et al. 



(120001 ). and calibrate them in terms of absolute (V — I) color by force- fitting them to the 
observed RGB sequences for M ilky Way globular clusters. This approach is described in 
detail in lHarris fc Harris! (120021 ) and we do not repeat it here. However, we emphasize that 
the metallicity scale we derive this way is an observationally calibrated one based on real 
globular clusters. The theoretical RGB tracks are used only to aid interpolation between 
the observed sequences for real clusters, and thus our results do not depend critically on the 
particular choice of stellar models. The one unavoidable and model-dependent assumption 
that underlies this method is that the mean age of the NGC 3377 stars is taken to be the same 
as for our local globular clusters (that is, about 12 Gy). If in fact they are systematically 
younger, then this method would slightly underestimate their Z— abundances since the RGB 
locus shifts blueward at lower age. However, as long as the age differences a re not severe, 
the age correction to the MDF is small; for example, in iRejkuba et al.l (120051 ) we note that 
a shift in mean age from 12 Gy to 8 Gy would produce only a 0.1— dex underestimat e in lo g 
Z. See also the RGB track comparisons for different ages in lHarris. Harris. &: Poold (119991 ). 
which demonstrate similar conclusions. 



In Figure H2] we show again the composite colo r-magnitude diagram, bu t now with the 
RGB tracks added. The solid lines are those of the IVandenBerg et al.l (120001 ) grid, whereas 
the two dashed lines at right are ones for Solar [Z = Z Q ) and ~ 3Z Q metallicities. It 
is clear, as we noted above, that the photometric detection limit for red stars set by the 
F606W exposures would eliminate any giant stars at Solar metallicity or above from our 
sample, even if they were present. What is not obvious from this Figure is whether or not 
this detection cutoff imposes any more basic limit on our interpretation of the MDF; that is, 
are many high— Z stars actually likely to be present in this galaxy? To gain a better reply to 
this question, we need to look at the details of the MDF itself as far as we can gauge them. 

Our derived MDF is shown in Figure [TBI where we divide the sample into half-magnitude 
bins by approximate luminosity Mboi- By doing this, we test for any systematic errors in the 
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interpolation procedure that might result from incorrect placement of the RGB model grid 
on the data; we also use it to reveal any broadening of the deduced MDF from photometric 
scatter toward the faint end of the data. Encouragingly, we see no systematic shift in the 
shape of the MDF: in all three bins, the peak occurs at log (Z/Z Q ) ~ —0.6, with a broad 
tail extending to lower metallicity and a much steeper ramp-down t o higher metallicity. As 



we found for the giant elliptical NGC 5128 ([Harris fc Harris! |2002| ). remarkably few stars 



are more metal-poor than log (Z/Z@) = —1.5. As a visual confirmation of this point, note 
directly from the CMD how few stars there are near the bluest RGB tracks, particularly 
at the high-luminosity end where the track color is most sensitive to metallicity. In other 
words, NGC 3377 - like the giant ellipticals, but unlike the smallest dwarfs and globular 
clusters - contains few stars indeed that are in the same regime as the "classic metal-poor" 
halo of the Milky Way. This comparison (see also the discussion below) indicates that NGC 
3377 did not form simply by the amalgamation of fully-formed dwarfs. 

Fig. [TBI explicitly shows the MDF with, and without, photometric completeness correc- 
tions, since these are important to our assessment of the MDF. Any extremely red stars 
falling to the right of the / =50% line in the CMD have been rejected from the sample, since 
in this region the completeness correction itself becomes dangerously large and the random 
and systematic errors of the photometry increase rapidly. For all stars brighter than this 
cutoff level, the completeness-corrected samples (the open histograms in Fig. [TB]) have been 
constructed by weighting each star individually as (1//) where / = fj ■ fy is the combined 
completeness fraction at its particular location in the CMD. For comparison, the unweighted 
MDF (based only on counting up all stars with / > 0.5) is shown in the hatched regions. 

Because the 50% completeness line in V runs nearly parallel to the RGB tracks in its 
region, it essentially has the effect of cutting off the MDF rather abruptly at log (Z/Z Q ) = 
—0.3. Redward of this point, we have no reliable knowledge of the MDF shape. But by 
the time we reach this photometrically-driven cutoff, the MDF is already well past its peak 
frequency and declining steeply. The evidence is thus strongly suggestive that we have seen 
the majority of the total MDF, and that only a small fraction of its stars are at Solar-type 
abundance or higher. By contrast, the giant NGC 5128 has an MDF extending well up past 



Solar abundance, even in the outer reaches of its halo ([Harris fc Harris! 120021 ; iRejkuba et al. 



20051) 



The fact that the V completeness cutoff is so close to the MDF peak, however, leaves 
lingering doubts that we have indeed seen the full shape of the metallicity distribution. Could 
there, in fact, be significant numbers of more metal-rich stars in the galaxy that perhaps 
belonged to a later star-forming episode and are still waiting to be identified? Such stars 
would be almost as bright in I as the other RGB stars, but would have been too red to have 
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been detected in V. To look into this a bit further, we have used the /—band measurements 
alone to isolate stars that were (a) detected and measured on the /—band image; (b) not 
detected on the V— band image; and (c) comparably bright in / to the other RGB stars. If 
there are large numbers of these, it would be suggestive that we have missed much of the 
actual RGB. 

The luminosity function for the stars measured only in /, and passing the same photo- 
metric culling described above, is shown as the dotted line in Figure dH This LF rises to a 
peak at the completeness limit / ~ 27.7 and then falls steeply. By contrast, the LF for the 
stars measured in both V and / (solid line in Fig. [HJ replotted from Fig. [TO]) has a broad 
peak determined by the V— band completeness limit. 

Very metal-rich RGB stars (those with log Z/Z Q > —0.3) would belong to the "/ only" 
LF and would lie in the approximate range from / ~ 26.6 (for the very top of the RGB at 
(V — I) ~ 3) down to the / completeness limit. We see that the LF for the stars in that range 
(dotted line in the Figure) stays well below the LF for the stars with both V, I (solid line). 
In the range 26.6 < I < 27.6, 20008 objects were measured only in / but not V. But this 
total should be considered as only an upper limit to the actual number of metal-rich RGB 
stars, because close visual inspection of the "/ only" candidates shows that a high fraction 
of them are small clumps of bad pixels, or faint, very small background galaxies. These 
passed the x an d &i selection steps, but did not match objects appearing at the same place 
on the V image and thus were rejected in the construction of the CMD. Taking these into 
account, we make a very rough estimate that there are fewer than ~ 10000 genuine RGB 
stars redder than the completeness cutoff in our CMD. By comparison, if we employ a simple 
chemical evolution model to predict the shape of the entire MDF and use it to extrapolate 
the observed MDF up past the photometric completeness cutoff (see next section), we find 
that 10% of the MDF should lie at Z > 0.5Z Q , corresponding to ~ 5700 'missed' stars in 
addition to our observed sample of 57000 over all radii. These very metal-rich stars would 
not change the MDF in a major way if (for example) they formed only the high-metallicity 
tail of the MDF that we already see. Unfortunately, any more definitive statement will 
have to await deeper photometry than is currently available. A clear answer would require 
doubling the V exposure time; that is, adding another 15 orbits to what we already have in 
hand. The incremental cost of HST time would thus be significant. 



5.2. A Chemical Evolution Model 



To step beyond the raw MDF i nto a simple phy s ical in terpretation, we use the basic 
chemical evolution model outlined in lHarris fc Harris! (120021 ) (hereafter HH02) that we ap- 
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plied to the NGC 5128 halo and bulge stars. Similar models have also been used for the 
halo of the Milky Way (|Prantzos 2003 ). and the globular cluster systems of large galaxies 



( jVandalfsen fc Harrisl 120041 ) . among other situations. This is the so-called "accreting box" 
model in which we assume that a region of primordial (Z ~ 0) gas turns itself into stars 
through a long succession of star-forming episodes, all the while that more gas is flowing 
into the region. The rate of gas infall is assumed to die away gradually and smoothly with 
time, so that in the late stages of the region's history, its chemical evolution asymptotically 
approaches the classic "closed box" or "simple" model. In a general sense, this model is a 
first-order description of what would be expected to happen during hierarchical merging of a 
large set of initial, pristine gas clouds, wherein star formation within the clouds is happening 
simultaneously with the merging into a bigger final galaxy. 

In HH02 we postulate that the gas infall rate starts at a rather high level and then 
dies away as an exponential decay with time. We also assume that at each star formation 
step, the same (small) fraction of ambient gas gets turned into stars. The abundance Z of 
the stars forming at any given moment is then determined by the Z(gas) left behind by 
the previous steps, in addition to the amount and composition of the new gas entering the 
region just before the next star formation step occurs. The sequence of timesteps can be 
numerically integrated to give the final distribution n(Z) of the stars (see HH02 for details). 
The model in its simplest form has three essential and unavoidable parameters: the effective 
yield y e ff of nucleosynthesis; the initial gas infall rate relative to the amount of gas initially 
present in the region; and the exponential decay time for the infall. 

An extremely effective way to exhibit the match between model and data is through 
the linear form of the MDF, as number of stars per unit heavy-element abundance Z/Z Q . 
The key results are shown in Figure [T5l where we include stars brighter than M^i ~ —2.5, 
i.e. ones in the upper two panels of Fig. [T3J As before, we include only stars at radii 
Rac > 2'. The raw (shaded) and completeness-corrected (unshaded) MDF are both shown. 
One particular accreting-box model solution is shown as the curved line. For Z > O.5Z , 
the observed MDF is cut off by photometric incompleteness, but as noted above, this cutoff 
is well past the MDF peak at Z ~ 0.15Z Q , and the MDF is seen to decline strongly up till 
the photometric cutoff. 

The model fit shown in Fig. [15] is for a triad of parameters that match the completeness- 
corrected data. The effective yield is y e ff = 0.0039 ~ O.26Z ; the exponential decay time 
for infall is 30 timesteps; and the total amount of gas infall is equal to 4 times the initial gas 
present in the region. The numerical integration is carried through several hundred timesteps 
until almost no gas remains. In each timestep we assume following HH02 that 5% of the 
ambient gas is converted to stars, and that the infalling gas has Z = 0. As we commented in 
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HH02, the absolute value of the timestep 5t is arbitrary in this model and must be supplied 
by other methods; however, through various arguments about the actual rate of hierarchical 
merging (see again HH02), it is plausibly near 5t ~ 30 Myr give or take factors of two. If 
so, the main period where infall is strongly influencing the MDF shape is thus the first 1 — 3 
Gyr. 

Each of the three critical parameters in the model controls three distinctly different 
features of the final MDF. First, the rapid rise in n(Z) from Z = to Z = 0.1 is set by the 
initial infall rate: the less the infall of pristine gas, the smaller the rise, until in the limit 
of zero infall, the n(Z) curve approaches the simple exponential decline of the closed-box 
model, n ~ e~ z l y . Second, the location of the peak, where the MDF turns over and begins 
to decline, is driven by the infall decay time: a longer decay time produces a higher— Z peak 
point since there is more gas to allow the enrichment process to continue driving upward. 
Finally, the long exponential tail to higher Z y where infall has essentially stopped and the 
only process left is the conversion of the remaining gas into stars (that is, a closed-box 
model), is shaped by the effective yield y e ff- 

We have restricted the present comparisons to the simplest set of parameters possible 
within the context of the model, and the three that are the least co-dependent. But additional 
model parameters are possible: for example, ones introduced in our previous analysis of NGC 
5128 included an initial, very short, time period T\ over which the gas infall is constant or 
rapidly ramps up; and also the assumed abundance Zj of the infalling gas. However, these 
prove to be much less critical to the solution than the three primary parameters defined 
above. The MDF data do not have enough resolution at extremely low Z to constrain n; 
and any plausibly small abundance for the infalling gas in the range Zi < 0.1Z Q can be used 
to generate adequate matches to the data (see HH02 and Vandalfsen & Harris 2004 for more 
extensive discussion). 

The value of the effective nucleosynthetic yield at about one-quarter Solar is a factor of 
two or three smaller than standard y— values that would hold in a closed-box model where all 
the gas remains held in situ and where the stellar IMF and nucleosynthesis rate are normal 
(e.g. iPagel &: Patchettl Il975l ). But if some fraction of the gas escape s in stellar and SN- 
driyen winds (leaky-box) then the effect is to lower y e ff in proportion (IBinney &: Merrifield 



19981 ). Formally, our results indicate that about half the gas present in the NGC 3377 region 
was ejected during the successive rounds of star formation, preventing the enrichment from 
proceeding upward past about one-half Solar levels. For comparison, in HH02 we found 
y e jf ~ 0.3Z G for the outer halo of NGC 5128, similar to the solution for NGC 3377. For 
the inner region of NGC 5128, much deeper in its potential well where virtually all the gas 
could be held, we found y e ff ~ 0.85Z Q . 
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For systems such as these, with broad MDFs and obviously extended star formation 
histories, much more complex models can be clearly brought to bear. Numerous examples 



of such models can be found in the literature (e.g. iNaab fc Ostriker 



2005 



2004 



Worthev et all 120051 ; IScannapieco et all 120051 ; IValle et al. 



2005 



2006 



Romano et al. 



Lafranchi fe Matteuci 



Oevlhoool . among many recent discussions). These have been developed all the way to 



full hierarchical-merging and chemodynamical-evolution codes (e .g. iKawata fc Gibson! 120031 ; 
Cole et al.ll2000l ; iBeasley et al.ll2003l ; ISomerville fc PrimackHl999 ). In so me cases very differ- 
ent parametric approaches are used; for example, IWorthey et all (120051 ) advocate the use of 
a closed-box model in which y e ff changes with time. In such a model, quite a strong decrease 
in y with increasing Z is necessary to reproduce the rising part of the MDF at low Z. In our 
view, the use of gas infall combined with y e ff ~ const is much more physically plausible, 
and has the strong advantage of employing a process (inflow of low-metallicity gas) which is 
almost certainly happening anyway at early times. 

It is notable, however, that the minimally simple set of parameters we use here already 
provides an accurate first-order fit to the MDF. The interpretation we offer for the validity 
of our first-order model is simply that in galaxies this large, we are looking at the combined 
results of so many star-forming events that much of the individual detail has been averaged 
out and is now no longer distinguishable. 



5.3. Gradients and Fine Structure? 



We can also use the MDF to search for any indication of a radial metallicity gradient. 
Our ACS/WFC field has stars with projected radial distances from the center of NGC 3377 
from 1'3 to 5'.7 (4.1 kpc to 17.7 kpc), a large enough range to be worth exploring for trends. 
The MDF, again in its linear form, is shown in Figure [TBI for an inner zone (1'.3— 2'0), a middle 
zone (2f0 — 3^5), and an outer zone (3'.5 — 5'7) . In each case, we show the same chemical 
evolution model applied to all three radial bins that we deduced for the total population in 
Fig. [T5] above, simply renormalized for the different total numbers of stars. 

All three panels are similar. We see no evidence for a radial gradient in metallicity, either 
in the peak position of the MDF or in its detailed shape. Although we deleted the innermost 
region R < 2' in the earlier analysis to minimize random scatter in the CMD, LF, and MDF, 
it has not led to any systematic bias in our understanding of the metallicity structure. This 
is, perhaps, not a surprising result since the effective radius of the galaxy is l'l and thus all 
the stars we sample are well outside it. In NGC 5128, we fo und that the MDF is virtually 
unchanged over a rad i al range from ~ 10 kpc out to 40 kpc ([Harris. Harris, fc Poolelll999 



Harris Harris! l2000i ; iRejkuba et al.ll2005l ); it is only within the bulge region at r < 8 kpc 
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(for NGC 5128, r e ff = 5.5 kpc) that significant changes began showing up and the MDF 
became systematically more metal-rich (HH02). 

Finally, we note that in either logarithmic or linear form, the MDF is very broad but 
contains a tantalizing hint of fine structure. The "notches" at (Z/Z Q ) ~ 0.2 and at ~0.3 are 
statistically significant at > 5cr in a nominal sense (the bin at the first gap has 1000 stars 
whereas the smooth model would predict 1500; and the second gap has 820 stars versus 
1120 expected from the model curve). Furthermore, the same gaps appear in all three of 
the radial bins plotted in Fig. [161 These small gaps are suggestive, perhaps, of distinct 
episodes of star formation that were incompletely averaged out over the whole galaxy. The 
extreme low-luminosity systems u Cen and Ca rina show MDFs that are more sharply divided 

among 



Koch et al. 2006; Sollima et al 



2005 



into a small number of distinct episodes (e.g. 
many others), whereas no such gaps appear in the much larger NGC 5128 (lHarris fc Harris 
20021 ; iRejkuba et al.ll2005l ). In this respect as in many others, NGC 3377 appears to occupy 
an intermediate position in the whole range of old, composite stellar systems: the smaller 
the system, the more likely it is that its MDF will be dominated by a smaller number of 
star-forming events and will take on a discontinuous morphology. 



5.4. Comparisons with Smaller and Larger Systems 

Adding this study to some of the previously cited papers, we can now make a first, 
admittedly rudimentary, attempt to put together a comparison of MDFs for the full range of 
normal E galaxies, from dwarfs to giants. An illustrative sample is shown in Figure [171 To 
represent the low-luminosity end of the normal E sequence, we use the "ty pical" Loca l Grou p 



dwarf NGC 147 (My = -15.6), with photometric (VJ) RGB data from lHan et all (119971 ). 



We also include two spatial regions from the giant NGC 5128, at My = —22.1: an outer-halo 



field at 40 kpc projected distance (IRejkuba et all 120051 ). and the 8-kiloparsec "inner halo 



or outer-bulge field analyzed by HH02. In all cases, exactly the same interpolation code and 
RGB model grid has been used to derive the MDF, so the results are internally homogeneous. 

The sequence in Fig. [T7] runs from the small potential well of the dwarf, through the 
intermediate-mass NGC 3377, to the outer regions of a giant, and finally to the deep inner 
regions of the giant. The change in the MDF clearly shows up as a progressively increasing 
spread in heavy-element abundance. We repeat that the upper tail for NGC 3377 is cut off 
by photometric incompleteness for Z > 0.QZ & and it is quite possible that a thin extension 
to higher Z is present. By contrast, the NGC 147 and NGC 5128 color-magnitude diagrams 
are more deeply sampled and less affected by incompleteness. In terms of physical trends, 
however, it is interesting to note that the peak of the MDF changes only modestly, from 
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~ 0.1Z Q in the dwarf up to ~ 0.3Z Q in the bulge of the giant. Instead, the most prominent 
effect shows up in the relative number of stars at high metallicity: in the dwarf, nothing is 
present above 0AZ Q , while the MDF for the giant extends far above Solar. This is the most 
dramatic effect of the depth of the potential well within which their stars were forming: in the 
deeper wells, almost all of the gas will be retained, permitting the enrichment to continue 
on through many generations and go up to high Z. A secondary feature worth noting is 
the steepness of the rise at low Z in the two smaller galaxies, showing the importance of 
infalling gas. This pattern is suggestive of an interpretation that the smaller galaxies took 
somewhat longer to assemble than the giants, during which time the inflow of pristine gas 
could continue longer and have a more important effect on the final MDF. 

6. Summary 

We have used deep HST/ACS images of a field in the intermediate-luminosity elliptical 
galaxy NGC 3377, a member of the Leo group, to derive the metallicity distribution function 
of its stars. More than 57000 stars were measured over a radial region extending from l'.Q 
to 5'7 from galaxy center. The (V, I) data reach deep enough to reveal almost 2 magnitudes 
of the red-giant branch and show the expected large spread in color that accompanies a 
high internal dispersion in metallicity. No significant presence of any young stellar popula- 
tion (stars younger than a few Gyr) can be seen in the color-magnitude diagram and the 
interpretation of the system is relatively simple. 

Using a finely spaced grid of RGB models calibrated against the Milky Way globular 
clusters, we derive the heavy-element abundance distribution n(Z) for the NGC 3377 stars. 
No apparent radial gradient in metallicity shows up; in all parts of the halo, the distribution 
peaks at log (Z/Z Q ) ~ —0.6 and is quite broad. Very few stars are, however, more metal- 
poor than log (Z/Z Q ) = —1.5, and the data suggest that few stars are more metal-rich than 
—0.3. It is possible that this galaxy did not reach Solar abundance during its enrichment 
history, at least in its halo. The overall shape of its MDF in all respects lies between the 
more metal-poor dwarf ellipticals on the one hand, and the very much more metal-rich giants 
on the other. 

WEH and GLHH thank the Natural Sciences and Engineering Research Council of 
Canada for financial support. We also thank an anonymous referee for constructive comments 
on the first version of the paper. 
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Table 1. Basic Parameters for NGC 3377 



Parameter 


Value 


Type 


E5 


a (J2000) 


10 fe 47 m 42 s4 


5 (J2000) 


13?59'08" 


v r (helio) 


665 km s _1 


Ay 


0.10 


(m - M) 


30.17 ±0.10 


V$ 


10.23 


Ml 


-19.9 



Fig. 1. — A reproduction of our measured ACS field in NGC 3377. Reproduced here is the 
combined F8UW (I band) image. At the ~ lf-Mpc distance of NGC 3377, the ACS/WFC 
field size corresponds to a width of 11 kpc. The inset box shows the location of the smaller 
field reproduced in the next Figure. Note that the center of the galaxy is well off the field 
to the lower right; the ACS field center is 3'8 away from the galaxy center. 
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Fig. 2. — A small section of the NGC 3377 ACS field. The image shown here is from the 
composite /—band exposure and is 36" (700 pixels) wide. The halo red-giant stars are the 
small starlike points scattered everywhere across the field; a few distant background galaxies 
are also visible. 
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Fig. 3. — Histograms of the nearest-neighbor distances (NND) for all detected objects on 
each of the two frames. We reject any objects with NND < 3 px from the subsequent 
analysis to minimize any effects of crowding (see text) 
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Fig. 4. — The color-magnitude array in (V, I) for our sample of halo stars in NGC 3377. All 
70,647 measured stars with high quality photometry are shown. The dashed lines show the 
50% completeness limits of our measurements, as derived from artificial-star tests (see text). 
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Fig. 5. — Photometric detection completeness functions, as evaluated by artificial-star exper- 
iments. The vertical axis plots the number of stars recovered by our DAOPHOT find/sdlstar 
measurement sequence divided by the number artificially inserted at the magnitude given 
along the horizontal axis. The solid symbols are the FS14W (/—band) stars and the open 
symbols are for FQOQW (V— band). The fitted curves drawn through each set of points are 
Pritchet functions with parameters as given in the text. 
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Fig. 6. — Color-magnitude diagrams subdivided by region. Here R is the distance from the 
center of NGC 3377 in arcminutes. 



- 30 - 




i i i i ! i j i i ! i i i i I i i i i I i i i i 

1 £ 3 4 5 6 

R cc (arcmin) 



f 

Fig. 7. — Plot of the local sky intensity level as a function of radius from the center of 
NGC 3377. Here the vertical scale is the local sky level around each star, in units of DU as 
returned by allstar. The vertical normalization is arbitrary. 
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Fig. 8. — Random measurement uncertainties of the photometry, as derived from the 
artificial-star tests. Each plotted point (in 0.02-magnitude bins) gives the root-mean-square 
scatter between the input and measured magnitude for all stars in the bin. The mean lines 
are exponential curves with parameters as given in the text. The two vertical dashed lines 
show the limiting magnitudes of the photometry, i.e. the magnitude levels at which the 
detection completeness is 50%. 
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Fig. 9. — Random and systematic errors of the photometry, as derived from artificial-star 
tests. Each graph shows the difference between the inserted and measured magnitudes of 
the fake stars in the sense (input minus measured), so that stars measured "too bright" have 
positive values in the plot. The solid line through each plot at 0.1-magnitude bins is the 
median Am within each bin. The vertical dashed line in each plot shows the "completeness 
limit" (the magnitude, from the previous figure, at which the detection completeness drops to 
50 percent). Systematic errors increase dramatically for stars fainter than the completeness 
limit. 
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Fig. 10. — Luminosity function for all measured stars in the NGC 3377 field. In the up- 
permost graph, the number of stars per unit magnitude is plotted against /; the faint-end 
steep decline is driven by strong photometric incompleteness. The middle panel shows our 
numerical estimate for the first derivative dn/dl of the LF, while the lower panel shows 
the numerical estimate for the second derivative d 2 n/dl 2 . In all panels the data have been 
convolved with a Gaussian smoothing kernel of a = 0.02 mag. 
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Fig. 11. — Left panel: Locations of stars in the brightest half-magnitude of the NGC 3377 
red-giant branch. The location of the gap between the two ACS/WFC chips can be seen 
through the middle, as well as small cutout masked regions around some background galaxies. 
Right panel: Locations of stars brighter than the RGB tip. 
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Fig. 12. — Color-magnitude diagram for the NGC 3377 red giants, with model red-giant 
tracks superimposed on the data points. All tracks are for ages of 12 Gy, but differ in metal- 
licity roughly in steps of A[Fe/H] ~ 0.1. These are the same set of tracks from VandenBerg 
et al. (2000) we have used for our previous study of the NGC 5128 halo, supplemented by 
two metal-rich tracks generated from old Milky Way star clusters (see Harris & Harris 2002), 
extending; from loe (Z/Z^) = —2.0 to +0.4. The Dlacement of the tracks assumes a fore- 
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Fig. 13. — Metallicity distribution function for the halo red giants, divided into three ap- 
proximate luminosity bins. This division tests for any systematic errors in the interpolation 
routine or the placement of the evolutionary tracks. The shaded regions show the MDF 
uncorrected for photometric completeness, while the higher unshaded regions show the full 
completeness-corrected MDF. Stars more metal-rich than log (Z/Z Q ) ~ —0.3 fall to the red 
of the 50% photometric completeness line (see previous Figure) and are not included in the 
sample. 
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Fig. 14. — Luminosity function for the NGC 3377 stars. The vertical axis here is rescaled 
from Figure 10 to show the number of stars per 0.05-mag bin plotted against / magnitude. 
The solid line shows the data for the stars in the previous figure, with measurements in both 
V and /. The dotted line shows stars measured only in / but not in V, thus ones that do 
not appear in the color-magnitude diagram. The vertical dashed line at / = 27.7 is the 50% 
photometric completeness level in I. Note that for red stars measured in both V and I the 
detection completeness is severely incomplete fainter than V ~ 26.6, thus the turnover in 
the solid curve is driven by the V filter cutoff. 
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Fig. 15. — Heavy-element abundance distribution, plotted in linear form. The number of 
stars per unit bin AZ = 0.02 is plotted against Z/Z Q . The solid curve is a chemical evolution 
model with exponentially declining infall, as discussed in the text. The assumed effective 
yield (the combination of stellar nucleosynthesis plus outflow; see text) is y e ff = 0.0039 ~ 
0.26Z m . 
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Fig. 16. — Heavy-element abundance distribution for three different radial ranges. The 
projected radial range in arcminutes is shown in each panel, from innermost to outermost. 
As in the previous Figure, the shaded region shows the raw number of stars in each Z— bin 
and the unshaded region the completeness-corrected bin totals. The solid curve in each case 
is the same chemical evolution model that was applied to the total population in the previous 
Figure. All three regions shown are outside the effective radius R e — l'.l of the galaxy. 
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Fig. 17. — Heavy-element abundance distributions for four different samples of stars. From 
top to bottom these are: (a) the Local Group dwarf elliptical NGC 147, with data from Han 
et al. (1997); (b) the Leo elliptical NGC 3377 (this study); (c) the outer halo of the giant 
elliptical NGC 5 128, from lReikuba et~aD tood ): and (d) the inner halo of NGC 5128, from 
Harris fc Harrisl koO$ ). 



